Discrete gap solitons in modulated waveguide arrays 
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We suggest a novel concept of diffraction management in waveguide arrays and 
predict the existence of discrete gap solitons that possess the properties of both 
conventional discrete and Bragg grating solitons. We demonstrate that both the 
soliton velocity and propagation direction can be controlled by varying the input 
light intensity. 
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Photonic structures with a periodic modulation of the 
refractive index (e.g. photonic crystals) can be used to 
precisely control propagation of optical pulses and beams. 
Wave localization is possible inside the band gaps of the 
linear spectrum, whereas dispersion and diffraction char- 
acteristics are strongly modified near the band edges. In 
particular, recent papersoatl reported on the fabrication 
of one-dimensional periodic arrays of identical optical 
waveguides where the effective diffraction coefficient can 
vanish or even become negative, being controlled by the 
input conditions and array parameters. 

In comparison with homogeneous media, efficiency of 
nonlinear processes can be greatly enhanced in prop- 
erly designed periodic structures. For waveguide arrays, 
where waves are primarily localized in weakly coupled 
waveguides, the effective diffraction can be greatly re- 
duced, further lowering the threshold for the beam self- 
focusing. It was predicted that diffractional spreading 
is suppressed for discrete solitonsux. 
possess many remarkable properties 
like their continuous counterparts, discrete solitons can 
propagate across an array at low powers, while at hish 
powers they become trapped by array discreteness, 
and this behavior is readily observed in experimental. It 
was also demonstrated that discrete solitons can be effi- 
ciently routed through two-dimensional networks of cou- 
pled waveguidescil 

In a number of recent studiesOJil, it was suggested 
that the properties of discrete solitons can be modified 
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FIG. 1: Schematic of a binary array of waveguides with 
alternating widths di and cfc and separation d s . 



by a periodic modulation of waveguides along the prop- 
agation direction. In this Letter, we suggest another, 
even more fundamental concept of the array engineering 
and consider a binary waveguide array composed of al- 
ternating "thick" and "thin" waveguides, as illustrated in 
Fig. [l]. In such a structure, the effective refractive index 
experiences additional transverse modulation and, there- 
fore, a "Rowland ghost gap" may appear in the linear 
spectruml!3. Formation of solitons in such gaps was ear- 
lier studied in the context of superstructure fiber Bragg 
gratingstJ, where the analysis is based on the coupled- 
mode equations. In this Letter, we demonstrate, for the 
first time to our knowledge, the existence of discrete gap 
solitons that display the properties of both conventional 
discrete and Bragg grating solitons and resemble nonlin- 
ear localized gap modes in diatomic latticesta. 

Propagation of waves in an array of weakly coupled 
single-made optical waveguides can be approximately de- 
scribedoEl by the discrete nonlinear Schrodinger equation 
(DNLS) for the normalized amplitude of the electric field 
E n localized at the waveguide with the index n, 



, dE n 
dz 



X n E n + (£„_i +E n+1 ) + ln \E n \ 2 E n = 0, (1) 



*URL: tittp: / /www.rsphysse . aim. edu . au/nonlinear 



where A„ characterizes the linear propagation constant of 
the mode guided by the n— th waveguide (which depends 
on its width), and 7„ is the effective nonlinear coeffi- 
cient. For the analysis of the structure shown in Fig. ^, 
it is convenient to introduce the notations: a n = E% n , 
b n = E 2n +i, A 2n +i = —p, and X 2n = 0, assuming the ap- 
propriate normalization of Eq. ([!]). In order to simplify 
the analysis, we neglect absorption and also consider the 
identical nonlinear coefficients, 7„ = 7. However, we 
have verified that the main conclusions of our study re- 
main valid if 7„ are weakly modulated. With no loss 
of generality, we have p > 0, so that a n and b n are the 
field amplitudes at the thick and thin waveguides, respec- 
tively. Dependence of the normalized detuning on the 
waveguide widths d± t 2 and the separation d s (see Fig. ^) 
can be found numerically. 

Our structure pan be compared to the waveguide ar- 
rays with defectsta, where localization at a thin waveg- 
uide {d,2 — 2.5 pm) embedded in an array of thicker 
waveguides (e?i = 4 pm) with d s — 5 pm was recently ob- 
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FIG. 2: Characteristic dependencies of the Bloch wave num- 
ber (Kb), group velocity (v g ), and the effective diffraction 
coefficient (D) on the propagation constant /3. Gray shadings 
mark the transmission bands. Normalized detuning between 
the thick and thin waveguides is p — 1 .5. 



served in experiments with AlGaAs arrays. Under such 
conditions, the normalized detuning between the modes 
of the thin and thick waveguides is about p ~ 1.5, the 
value we use in numerical simulations presented belpw. 
Defects with other parameters can also be fabricatedtZl. 

First we analyze the properties of linear Bloch waves 
in such a periodic binary structure. The Bloch waves are 
characterized by the wave number Kb, 



a„ = Ae if3z+iK » n , b n = Be i0z+iK >> n , 



(2) 



where [3 is the mode propagation constant. We substitute 
Eq. (j^) into the linearized equation ([j]) (with 7 = 0), 
and obtain the linear dispersion relation that couples the 
propagation constant and Bloch wave number, and yields 
a relation between the amplitudes at the thin and thick 
waveguides, 



K b = cos- 1 (n/2), A[3 = Be- iK "l 2 y/2- 



(3) 



where 77 = 2 — (3(j3 + p). It follows that the transmission 
bands, corresponding to real Kb, appear when /?_ < j3 < 
-p or < (3 < f3+, where (3+ = -(p/2) ± ^J(p/2) 2 + 4. 
On the other hand—an additional gap, also called the 
Rowland ghost gaplj, appears for — p < (3 < 0, and it 
increases for a larger difference between the widths of 
the neighboring waveguides. A characteristic dispersion 
relation and the corresponding band-gap structure are 
presented in Fig. ||(a). 



Using the dispersion relation (||), we calculate the 
group velocity, v g = —2(d(3/dK b ), and the effective 
diffraction coefficient, D = —4(d 2 f3/dK 2 ), of the binary 
array. Characteristic dependencies of these parameters 
on the propagation constant (3 are shown in Figs. ^|(b,c). 
In both transmission bands, there exist regions with ef- 
fective normal and anomalous diffraction separated by a 
zero-diffraction point. 

When the medium possesses a self-focusing nonlinear- 
ity, bright discrete solitons can form provided the effec- 
tive diffraction coefficient is positive. We consider exci- 
tation of such solitons by an input Gaussian beam, 

E n (z = 0) = Ccxp [-(n - n ) 2 /v 2 + w(n - n )] , (4) 

where no is the position of the beam center, v is the beam 
width, and k characterizes the inclination angle. Such an 
input beam can be presented as a superposition of two 
modulated linear eigenmodes, whose eigenvalues [3\^ 2 cor- 
respond to the Bloch wave number Kb — cos -1 [cos(2k)], 

E 2n (Q) = (^1 + A 2 ) exp [-(2n - n ) 2 /is 2 + iK b n] , 

E 2n +i(0) = (B 1 + B 2 ) exp [-(2n + 1 - n Q ) 2 /v 2 + iK b n] 

where the amplitudes Aj and Bj (j = 1, 2) satisfy Eq. (0) 
at (3 = [3j (we choose [3\ > f3 2 , with no loss of generality) . 
For a wide input beam (v > 4), a ratio of the powers of 
two excited Bloch modes can be written in a simple form, 



P 2 



(5) 



where 6 



sj\ + p/fk cosk /cos(K b /2). The two 

modes always have opposite group velocities and diffrac- 
tion coefficients (see Fig. ||) . 

It is useful to compare the results with a homogeneous 
array, when p = 0. In this case, within the validity of the 
DNLS model, there exists only one parameter domain 
where diffraction is positive, n a»d discrete solitons can 
exist, for self- focusing mediald'BB. On the other hand, 
in a modulated array diffraction is modified close to the 
Rowland ghost gap, and there appear two regions with 
positive diffraction, where solitons can exist. Thus, the 
binary array structure allows the band gap engineering 
in a broad parameter region. 

Let as first discuss the excitation of discrete solitons by 
a beam incident at a normal angle, so that k = Kb = 0, 
and /?i t 2 = /?+,-• In this case, p < 1, i.e. the first 
(j = 1) mode is always dominant (for example, p ~ 0.03 
for p = 1.5). As is demonstrated in Fig. g, at [3 — (3+ 
the effective diffraction coefficient is positive and, there- 
fore, discrete solitons can be formed in a self-focusing 
medium. Properties of such solitons are similar to those 
in homogeneous arrays. 

On the other hand, discrete gap solitons may appear 
due to localization withing the linear gap (at — p < [3 < 
0), which appears due to a peridic modulation of waveg- 
uide widths as discussed above. In order to excite sta- 
tionary gap solitons in a self-focusing medium, one might 
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FIG. 3: Excitation of discrete gap solitons in a self-focusing 
medium (7 = +1) by input beams with k = 0.6-7T, no = 1, v = 
4, and normalized peak intensities (a) |C| 2 = 0.3, (b) 0.75, 
(c) 0.9, (d) 5. The parameters correspond to Fig. 



attempt to use an input beam with k = tt/2, so that 
Kb = 7T, Pi = 0, fa — ~P- I n this case, two Bloch modes 
have equal powers (p = 1), however one of them expe- 
riences self-focusing and another one — self-defocusing, 
and in this situation an efficient generation of gap solitons 
is not possible. 

The optimum conditions for generating discrete gap 
solitons can be realized when tt/2 < k < it, and -D(/?2) > 
0. Indeed, under such conditions p > 1, i.e. the sec- 
ond mode, which experiences self-focusing, is dominant. 
Moreover, the Bloch wave envelopes have opposite group 
velocities, so that they move apart in the opposite direc- 
tions. In Fig. ||(a), we show that two beams are indeed 
generated at the input. The beam which moves to the 
right is localized at odd (i.e. thin) waveguides, and it 
transforms into a gap soliton. On the contrary, the other 
beam moves to the left, and it experiences self-defocusing 
and broadens. 

Stationary gap solitons can be generated at higher in- 
put intensities. Under such conditions, two Bloch modes 
are initially trapped together, resulting in a periodic 
beating that can strongly affect the soliton formation pro- 
cess. If the emerging soliton has a small velocity, it is not 
able to overcome the Peierls-Nabarro potential of the pe- 
riodic structure, and it becomes trappedEl. In this case, 
the gap soliton can remain in the center, as is shown 
in Fig. 0(b). However, if the soliton acquires a higher 



velocity, it can still move across the array, as is shown 
in Fig. y(c). Somewhat surprisingly, in the latter case 
the soliton moves to the left, i.e. in the direction oppo- 
site to the propagation direction of the input beam. If 
the intensity is increased even further, then the gap soli- 
tons become oscillatory-unstable, similar to the solitons 
in fiber Bragg gratingslia [see Fig. |](d)]. 

In conclusion, we have studied the diffraction proper- 
ties and nonlinear wave propagation in binary waveguide 
arrays with alternating waveguide widths. We have pre- 
dicted the existence of discrete gap solitons and demon- 
strated their intriguing dynamics controlled by varying 
the input intensity. 

The authors thank J. S. Aitchison, R. Morandotti, and 
Y. Silberberg for useful discussions. This work was sup- 
ported by the Australian Research Council. 
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